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Abstract 

Differential cross sections for jet photoproduction in the reaction ep — > e jet X 
have been measured with the ZEUS detector at HERA using 82.2 pbl of in- 
tegrated luminosity. Inclusive jet cross sections are presented as a function of 
the jet transverse energy, for jets with > 17 GeV and pseudorapidity 
— 1 < rf ct < 2.5, in the 7p centre-of- mass-energy range 142 < W 1P < 293 GeV. 
Scaled jet invariant cross sections are presented as a function of the dimensionless 
variable x T = 2E^/W 1P for (W 1P ) = 180 and 255 GeV. Next-to-leading-order 
QCD calculations give a good description of the measured differential cross sec- 
tions in both magnitude and shape. The ratio of scaled jet invariant cross sections 
at the two (W w ) values shows clear non-scaling behaviour. A value for the strong 
coupling constant of a s (M z ) = 0.1224 ± 0.0001 (stat.) ±g:8oM ( ex P-) -omit ( th 
has been extracted from a QCD analysis of the measured da/dE^ 1 . The varia- 
tion of a s with E^ 1 is in good agreement with the running of a s as predicted by 
QCD. 
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1 Introduction 



Jet production provides a testing ground for the theory of the strong interaction between 
quarks and gluons, namely quantum chromodynamics (QCD). This letter concentrates 
on one aspect of jet production, namely, the comparison of jet cross sections for the 
same reaction at different centre-of-mass energies. This highlights the effects of scaling 
violations, while a QCD analysis of jet-production rates allows the measurement of the 
strong coupling constant, a s . 

The parton model predicts a jet cross section that scales with the centre-of-mass energy. In 
this case, the scaled jet invariant cross section, (E 1 ^ YE* et d?a / dp i xdp i ydp i ^' , as a function 
of the dimensionless variable xt = 2E 3 ^ t /W, should be independent of W, where W is 
the centre-of-mass energy, E^ et is the jet energy, E^ is the jet transverse energy and 
(Px iPy \^z") are the components of the jet momentum. Thus, the ratio of scaled jet 
invariant cross sections for different centre-of-mass energies will be unity for all xt- On the 
other hand, QCD predicts that jet cross sections should exhibit a non-scaling behaviour, 
due both to the evolution of the structure functions of the colliding hadrons and to the 
running of a s . Scaling violations have been observed in the ratio of the scaled jet invariant 
cross sections as a function of xt in pp collisions at centre-of-mass energies of either 546 
or 630 and 1800 GeV [1]. 

At HERA, similar tests can be made in the photoproduction of jets. Two types of QCD 
processes contribute to jet production in jp interactions at 0(aa s ) [2,3]: either the photon 
interacts directly with a parton in the proton (the direct process) or the photon acts as 
a source of partons, one of which interacts with a parton in the proton (the resolved 
process). Violations of scaling should be observed both in resolved and direct processes. 
Furthermore, measurements of high-E 3 ^ 1 jet cross sections in 'jp interactions over a wide 
range of E^ allow a determination of a s (Mz) as well as its energy-scale dependence. 

This letter presents a measurement of the inclusive jet cross section in 7p interactions as 
a function of Ej, in the ^p centre-of-mass-energy range 142 < W 1P < 293 GeV for jets 
with pseudorapidity — 1 < rf ct < 2.5. Scaled jet invariant cross sections are also presented 
as a function of xt for (W JP ) = 180 and 255 GeV in the region —2 < rf^ < 0, where rf^ 
is the jet pseudorapidity in the 7p centre-of-mass frame. 

2 Experimental conditions 

The data were collected during the running period 1998-2000, when HERA operated with 
protons of energy E p = 920 GeV and electrons or positrons of energy E e = 27.5 GeV, and 
correspond to an integrated luminosity of 82.2 ± 1.9 pbl. A detailed description of the 
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ZEUS detector can be found elsewhere [4,5]. A brief outline of the components that are 
most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [6], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift-chamber layers, organized in nine superlayers covering the 
polar-angleQ region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as cr(p T )/p T = 0.0058^ © 0.0065 © 0.0014/p T , with p T in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the 
energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [7] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section (EMC) and either one (in RCAL) or two (in 
BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter 
is called a cell. Under test-beam conditions, the CAL single-particle relative energy 
resolutions were a{E)/E = 0.18/ y/E (GeV) for electrons and o{E)/E = 0.35/ y/E (GeV) 
for hadrons. 

The luminosity was measured from the rate of the bremsstrahlung process e + p —>■ e +/ -fp. 
The resulting small-angle energetic photons were measured by the luminosity monitor [8], 
a lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3 Data selection and jet search 

A three-level trigger system was used to select events online [5,9]. At the first level, events 
were triggered by a coincidence of a regional or transverse energy sum in the CAL and at 
least one track from the interaction point measured in the CTD. At the second level, a 
total transverse energy of at least 8 GeV, excluding the energy in the eight CAL towers 
immediately surrounding the forward beampipe, was required, and cuts on CAL energies 
and timing were used to suppress events caused by interactions between the proton beam 
and residual gas in the beampipe. At the third level, a jet algorithm was applied to 
the CAL cells and jets were reconstructed using the energies and positions of these cells. 
Events with at least one jet with > 10 GeV and rj < 2.5 were accepted. 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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Events from collisions between quasi-real photons and protons were selected offline using 
similar criteria to those reported in a previous publication [10]. The main steps are briefly 
discussed here. After requiring a reconstructed event vertex consistent with the nominal 
interaction position and cuts based on the tracking information, the contamination from 
beam-gas interactions, cosmic-ray showers and beam-halo muons was negligible. Charged 
current deep inelastic scattering (DIS) events were rejected by requiring the total missing 
transverse momentum, p™ lss , to be small compared to the total transverse energy, E^, 
i.e. p™ lss / \fl^ < 2 \/GeV. Any neutral current (NC) DIS events with an identified 
scattered-positron or electron candidate in the CAL [11] were removed from the sample 
using the method described previously [12]. The remaining background from NC DIS 
events was estimated by Monte Carlo (MC) techniques to be below 0.3% and was ne- 
glected. The selected sample consisted of events from ep interactions with Q 2 < 1 GeV 2 
and a median Q 2 pa 10~ 3 GeV 2 , where Q 2 is the virtuality of the exchanged photon. The 
events were restricted to 7p centre-of-mass energies in the range 142 < W 1P < 293 GeV, 
as described in Section |j. 

The longitudinally invariant fcy cluster algorithm [13] was used in the inclusive mode [14] 
to reconstruct jets in the hadronic final state from the energy deposits in the CAL cells. 
The jet search was performed in the pseudorapidity-azimuth (r)—(p) plane of the laboratory 
frame. The jet variables were defined according to the Snowmass convention [15]. The 
jets reconstructed from the CAL cell energies are called calorimetric jets and the variables 
associated with them are denoted by E^ l cal , r?^ and A total of 197155 events with 
at least one jet satisfying Ej?^ > 13 GeV and —1 < rf^ < 2.5 were selected. 

4 Monte Carlo simulation 

The MC programs PYTHIA 6.1 [16] and HERWIG 5.9 [17] were used to generate resolved 
and direct photoproduction events. In both generators, the partonic processes are simu- 
lated using leading-order matrix elements, with the inclusion of initial- and final-state par- 
ton showers. Fragmentation into hadrons was performed using the Lund string model [18] 
as implemented in JETSET [19] in the case of PYTHIA, and a cluster model [20] in the 
case of HERWIG. The generated events were used for calculating energy and acceptance 
corrections. The corrections provided by PYTHIA were used as default values and those 
given by HERWIG were used to estimate the systematic uncertainties coming from the 
treatment of the parton shower and hadronisation. Samples of PYTHIA including multi- 
parton interactions [21] with a minimum transverse momentum for the secondary scatter 
of 1 GeV [22] were used to study the effects of a possible "underlying event" . 

All generated events were passed through the ZEUS detector- and trigger-simulation pro- 
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grams based on GEANT 3.13 [23]. They were reconstructed and analysed by the same 
program chain as the data. The jet search was performed using the energy measured in 
the CAL cells in the same way as for the data. The same jet algorithm was also applied 
to the final-state particles; the jets found in this way are referred to as hadronic jets. 

5 Fixed-order QCD calculations 

The QCD calculations, at both leading order (LO) and next-to-leading order (NLO), 
used in this analysis are based on the program by Klasen, Kleinwort and Kramer [24]. 
The calculations use the phase-space-slicing method [25] with an invariant-mass cut to 
isolate the singular regions of the phase space. The number of flavours was set to five; the 
renormalisation, fj,R, and factorisation scales, /i^, were set to /ir = /i^ = [i = E^ ', a s was 
calculated at two loops using A^ = 220 MeV, which corresponds to a s {M z ) = 0.1175. 
The MRST99 [26] parameterisations of the parton distribution functions (PDFs) of the 
proton and the GRV [27] sets for the photon were used as defaults for the comparisons 
with the measured cross sections. 

Since the measurements refer to jets of hadrons, whereas the QCD calculations refer to 
partons, the predictions were corrected to the hadron level using the MC models. The 
multiplicative correction factor, Chad; defined as the ratio of the cross section for jets of 
hadrons over that for jets of partons, was estimated with the PYTHIA and HERWIG 
programs. The values of Chad obtained with PYTHIA were taken as the defaults; the 
predictions from the two models were in good agreement. The values of Chad differed 
from unity by less than 2.5%. 

6 Energy and acceptance corrections 

The comparison of the reconstructed jet variables for the hadronic and the calorimetric jets 
in simulated events showed that no correction was needed for "rf eXl and (fj ot (r]i ct ~ rf^i and 
yjjct ^ ip 3 ^)- However, the transverse energy of the calorimetric jet was an underestimate 
of the corresponding hadronic jet energy by an average of ~ 15%, with an r.m.s. of 
~ 10%. This underestimation was mainly due to the energy lost by the particles in the 
inactive material in front of the CAL. The transverse-energy corrections to calorimetric 
jets, as a function of rj^ and -Ej^ cal and averaged over were determined using the 
MC events. Henceforth, jet variables without subscript refer to the corrected values. 
After these corrections to the jet transverse energy, events with at least one jet satisfying 
E]£ > 17 GeV and —1 < rf ct < 2.5 were retained. 
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The 7p centre-of-mass energy is given by W 1P = y^sy, where y is the inelasticity variable 
and \/i is the ep centre-of-mass energy, s = 4E e E p . The inelasticity variable was recon- 
structed using the method of Jacquet-Blondel [28], yje — (E — pz)/2E e , where E is the 
total CAL energy and pz is the Z component of the energy measured in the CAL cells. 
The value of y was systematically underestimated by ~ 20% with an r.m.s. of ~ 10%. 
This effect, which was due to energy lost in the inactive material in front of the CAL and 
to particles lost in the rear beampipe, was satisfactorily reproduced by the MC simulation 
of the detector. The MC event samples were therefore used to correct for this underesti- 
mation [29] and obtain y cor . Events were required to have 142 < W 1P < 293 GeV, where 

The variable xt was reconstructed using the corrected values of Ej, and W 7P . Its resolu- 
tion was ~ 12%. The variable rf^ was computed by boosting ??> et to the jp centre-of-mass 
frame using the formula r]^ p = rf et — ]n(2E p /Wr fp ). The comparison of rf^ for the hadronic 
and the calorimetric jets in simulated events showed a good correlation, so that no cor- 
rection was needed. The resolution on rff p was ~ 0.08. 

The PYTHIA MC event samples of resolved and direct processes were used to compute the 
acceptance corrections to the jet distributions. These correction factors took into account 
the efficiency of the trigger, the selection criteria and the purity and efficiency of the jet 
reconstruction. The contributions from direct and resolved processes in the MC models 
were added according to a fit to the uncorrected data distribution of the energy deposited 
in the RCAL. A reasonable description of the E^, 77 jet , W 1P , r/^ and xt distributions in 
the data was provided by both PYTHIA and HERWIG. The differential inclusive jet cross 
sections were obtained by applying bin-by-bin corrections to the measured distributions. 
These correction factors differed from unity by typically less than 10%. 

7 Experimental uncertainties 

A detailed study of the experimental systematic uncertainties of the cross-section mea- 
surements included the following sources: 

• the effect of the presence of a possible underlying event was estimated by using the 
samples of PYTHIA including multiparton interactions to evaluate the correction fac- 
tors. This effect was typically 5% and increased to ~ 10% in the high-x^ tail of the 
scaled jet invariant cross sections; 

• the effect of the treatment of the parton shower and hadronisation was estimated by 
using the HERWIG generator to evaluate the correction factors. The uncertainty in 
the cross sections was typically 2%; 
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• the effect of the uncertainty on the modelling of the Q 2 spectrum of resolved pro- 
cesses in the MC was estimated by using the different approximations implemented in 
PYTHIA and HERWIG. The uncertainty in the cross sections was below 2%; 

• the effect of the uncertainty on W 1P was estimated by varying ?/jb by ±1% in simulated 
events. The uncertainty in the cross sections was below 1% at low E^ 1 , increasing to 
~ 3% at high E^; 

• the effect of the uncertainty on the parameterisations of the proton and photon PDFs 
was estimated by using alternative sets of PDFs in the MC simulation to calculate 
the correction factors. The variation of the cross sections was smaller than 1% in each 
case. 

The uncertainty on the simulation of the trigger was negligible. All the above systematic 
uncertainties were added in quadrature, giving a total systematic uncertainty in the cross 
sections of 5% at low E^ , increasing to ~ 10% at high E^ . The absolute energy scale 
of the calorimetric jets in simulated events was varied by its uncertainty of ±1% [10,30]. 
The effect of this variation on the inclusive jet cross sections was typically =f5% at low 
E^ increasing to =f10% at high E^. This uncertainty is highly correlated between 
measurements in different bins. The uncertainty in the luminosity determination of 2.25% 
was not included. 

8 Uncertainties on the theoretical predictions 

The following uncertainties were considered: 

• the uncertainty on the NLO calculations due to higher-order terms was estimated by 
varying /i between E^ /2 and 2E%\ It was less than 10% and mainly affected the 
normalisation. In the ratio of the scaled jet invariant cross sections, it was less than 
2.5%; 

• the uncertainty on the NLO calculations due to the uncertainties on the photon PDFs 
was estimated by using an alternative set of parameterisations, AFG-HO [31]. The 
effect was below 5% for the cross sections and 2% for the ratio; 

• the uncertainty on the NLO calculations due to the statistical and correlated system- 
atic experimental uncertainties of each data set used in the determination of the proton 
PDFs was calculated, making use of the results of an analysis [32] that provided the 
covariance matrix of the fitted PDF parameters and the derivatives as a function of 
Bjorken x and /if,. The resulting uncertainty in the cross sections was 1% at low Ei? 
and increased to 5% at high E^ . The uncertainty in the ratio of the scaled jet invari- 
ant cross sections was below 0.3%. To estimate the uncertainties on the cross sections 
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due to the theoretical uncertainties affecting the extraction of the proton PDFs, the 
calculation of all the differential cross sections was repeated using a number of differ- 
ent parameterisations obtained under different theoretical assumptions in the DGLAP 
fit [32]. This uncertainty was below 3% for the cross sections and negligible for the 
ratio; 

• the uncertainty on the NLO calculations due to that on a s [Mz) was estimated by 
varying a s (Mz) within its uncertainty [33] and, simultaneously, by repeating the cal- 
culations using two additional sets of proton PDFs, MRST99|T an d MRST99||, de- 
termined assuming a s (M z ) = 0.1225 and 0.1125, respectively. The difference between 
the calculations using these sets and MRST99 was scaled by 60% to reflect the current 
uncertainty on the world average of a s [33]. The resulting uncertainty in the cross 
sections was ~ 8% at low E 3 ^ decreasing to ~ 2% at high E 3 ^ 1 . The uncertainty in 
the ratio of the scaled jet invariant cross sections was below 4%; 

• the difference in the hadronisation corrections as predicted by PYTHIA and HERWIG 
resulted in an uncertainty smaller than 2.5%. 

All the above theoretical uncertainties were added in quadrature. 

9 Results 

9.1 Inclusive jet differential cross sections 

Using the selected data, inclusive jet differential cross sections were measured for 142 < 
W 1P < 293 GeV. The cross sections were determined for jets with E 3 ^ > 17 GeV and 
— 1 < r] 3Ct < 2.5. There were 113 843 events, containing 145 797 jets, in this kinematic 
region. 

The cross-section da/dE 3 ^, measured in the E 3 ^ range between 17 and 95 GeV, is pre- 
sented in Fig. p] and Table [l]. The data points are located at the weighted mean of each E 3 ^ 
bin. The measured da/dE 3 ^ 1 falls by over five orders of magnitude in this E 3 ^ range. Fig- 
ure^ and Table ^| show the scaled jet invariant cross-section, (E 3 ^ t ) 4: (E 3et d 3 a/dp 3 ^'dp l Ydp 3 ^ t ) 
averaged over the range —2 < r& < 0, as a function of xt for (W w ) values of 180 and 
255 GeV; the (W w ) values were chosen as the centres of the intervals 169-191 GeV and 
240-270 GeV. The measurements were restricted to the same range in to have the 
same acceptance for the two (W 7P ) intervals. 

Fixed-order QCD calculations are compared to the data in Fig. |l|. The LO QCD calcu- 
lation underestimates the measured cross section by about 50% for E 3 ^ < 45 GeV. The 
calculation that includes NLO corrections gives a good description of the data within 
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the experimental and theoretical uncertainties over the complete E 3 ^ range studied. In 
particular, no significant deviation is observed in the highest E 3 ^ 1 region. The NLO calcu- 
lations also give a good description of the scaled jet invariant cross sections as a function 
of XT-, as shown in Fig. [| 



9.2 Test of scaling 

To test the scaling hypothesis, the ratio of the scaled jet invariant cross sections as a 
function of xt was measured for the two chosen values of (W 1P ), after correcting for the 
difference in the photon flux [34] between these intervals. Figure ^| and Table |2| show the 
measured ratio as a function of xt- It shows a clear deviation from unity, in agreement 
with the NLO QCD predictions, which include the running of a s and the evolution of 
the PDFs with the scale. This constitutes the first observation of scaling violations in 7^ 
interactions. 

The ratio of the scaled jet invariant cross sections can be used to test QCD more precisely 
than is possible with the individual cross sections, since the experimental and theoretical 
uncertainties partially cancel. In particular, the experimental uncertainty in the absolute 
energy scale of the jets cancels almost completely in the ratio. The theoretical uncertainty 
on the predictions of the scaled jet invariant cross section was 13%, whereas that on 
the ratio was reduced to 2 — 5%. The NLO QCD prediction is in agreement with the 
data within the improved experimental (below 12%) and theoretical uncertainties. This 
agreement shows that the energy-scale dependence predicted by QCD is in accord with 
the measured dependence. 



9.3 Determination of cx s (Mz) 

The measured cross-section da/dE^ as a function of E 3 ^ was used to determine a s (Mz) 
using the method presented previously [35]. The NLO QCD calculations were performed 
using the three MRST99 sets of proton PDFs, central, MRST99|| and MRST99TT; the 
value of a s {Mz) used in each partonic cross-section calculation was that associated with 
the corresponding set of PDFs. The a s (M z ) dependence of the predicted da/dE 3 ^ in 
each bin i of E 3 ^ was parameterised according to 



da/dE2?(a.(M z )) . = C\a s [M z ) + C l 2 a 2 s {M z 



where C\ and C\ are constants, by using the NLO QCD calculations corrected for hadro- 
nisation effects. Finally, a value of a s [Mz) was determined in each bin of the measured 
cross section as well as from all the data points by a \ 2 fit- 
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The uncertainties on the extracted values of a s (Mz) due to the experimental systematic 
uncertainties were evaluated by repeating the analysis for each systematic check presented 
in Section [?[ The largest contribution to the experimental uncertainty comes from the jet 
energy scale and amounts to ±1.5% on a s (Mz)- The theoretical uncertainties were eval- 
uated as described in Section |[ The largest contribution was the theoretical uncertainty 
on a s (Mz) arising from terms beyond NLO, which was ^.33%- The change of a s (Mz) 
due to the uncertainties on the photon PDFs and on the hadronisation corrections were 
+0.7% and +0.8%, respectively. The uncertainty on a s {Mz) due to the uncertainties on 
the proton PDFs was ±0.9%. The total theoretical uncertainty on a s [Mz) was obtained 
by adding these uncertainties in quadrature. 

The values of a s (Mz) as determined from the measured dajdE 3 ^ in each region of E 3 ^ 
are shown in Fig. |4]a) and Table |3|. By combining all the E 3 ^ regions, the value of a s [Mz) 
obtained is 

a s (M z ) = 0.1224 ±0.0001 (stat.) ±g;g™ ( ex P-) -omit (th.). 

This value of a s [Mz) is consistent with the current world average [33] of 0.1183 ± 0.0027 
as well as with recent determinations from jet production in NC DIS at HERA [36, 35] 
and pp collisions at Tevatron [37]. It has a precision comparable to the values obtained 
from e + e~ interactions [33]. 



9.4 Energy-scale dependence of ot s 

The QCD prediction for the energy-scale dependence of the strong coupling constant 
was tested by determining a s from the measured do/dE^ at different E 3 ^ values. The 
method employed was the same as described above, but parameterising the a s dependence 
of dajdE 3 ^ in terms of a s ((E 3 ^ t )) instead of a s {Mz), where (Ej?) is the weighted mean 
of Ej? in each bin. The measured a s (E 3 ^ t ) values are shown in Fig. f|b) and Table |3|. The 
results are in good agreement with the predicted running of the strong coupling constant 
over a large range in E 3 ^. 

The energy-scale dependence of the measured a s (E 3 ^ t ) was quantified by fitting the results 
using the functional form predicted by the renormalisation group equation. Perturbative 
QCD predicts that a~\E^) varies approximately linearly with laE 3 ^ . At two loops, the 
energy-scale dependence of a~ 1 (E 3 ^ t ) is given by 

/? 1 ln(21n(4 et /A)) l" 1 
P 2 ln(4 Gt /A) J ' U 

where (3q — 11 — In/, /3i = 51 — and rif is the number of active flavours. Thus, the 
slope of a~ l {E 3 r^) gives I3q/2tx. A x 2 At to the extracted a~ 1 (E 3 ^ t ) values was performed 



«; 1 (4 et ) = SM4 et /A)- 
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to determine (3q using the functional form given by Eq. ([!]), leaving (3q and A as free 
parameters; (3\ was set to (19/?o — 107)/2. Figure |c) shows the measured a~ 1 {E i ^ t ) 
as a function of lni^ together with the results of the fit. Although the value of A, 
0.535 ± 0.073 (stat.) ±g;^jj (exp.) +°;^ (th.) GeV, is not well constrained in the fit, 
the value of a s (Mz) obtained by extrapolation from the results of the fit is more precise, 
a s (M z ) = 0.1188±0.0009 (stat.) (exp.) ±8:Sot? ( th 0- Thi s determination of a s (M z ) 

is consistent with that of Section |9.3|, in which the running of a s as predicted by QCD 



was assumed. The extracted value of (3q is 

fa = 8.53 ±0.22 (stat.) (exp.) ^ (th.). 

This value is consistent with the prediction of perturbative QCD for the relevant number 
of active flavours in the E 3 ^ 1 region considered, (3q = 7.67 for rif = 5. 



10 Summary and conclusions 

Measurements of differential cross sections for inclusive jet photoproduction have been 
made in ep collisions at a centre-of-mass energy of 318 GeV using 82.2 pbl of data 
collected with the ZEUS detector at HERA. The cross sections refer to jets identified 
with the longitudinally invariant kx cluster algorithm in the inclusive mode and selected 
with > 17 GeV and — 1 < rf ct < 2.5. The measurements were made in the kinematic 
region defined by Q 2 < 1 GeV 2 and 142 < W w < 293 GeV. 

The inclusive jet cross section was measured as a function of E^ in the range between 
17 and 95 GeV. The scaled jet invariant cross sections, averaged over —2 < < 0, were 
measured as a function of the dimensionless variable xt for (W 1P ) = 180 and 255 GeV. 
The NLO QCD calculations give a good description of the shape and magnitude of the 
measured cross sections. No significant deviation with respect to QCD was observed up 
to the highest scale studied. The ratio of scaled jet invariant cross sections at two values 
of (W 1P ) represents the first observation of scaling violations in 7p interactions. 

A QCD analysis of the measured da/dE^ 1 yields a value of the strong coupling constant 
of 

a s (M z ) = 0.1224 ±0.0001 (stat.) t° S ( ex P-) -aoo42 (th-), 

which is in agreement with the current world average and constitutes the first determina- 
tion of a s (Mz) from jet production in 7p interactions. The value of a s as a function of 
Ej, is in good agreement, over a wide range of E^ , with the running of a s as predicted 
by QCD. 
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(E8?) (GeV) 


da/dE 3 ^ ± stat. ± syst. (pb) 


syst. -scale (pb) 


18.6 


290 ± 1 ± 14 


(+11,-12) 


22.7 


97 ± 1 ±4 


(±4, -4) 


26.7 


37.6 ±0.3 ±2.6 


(+1.6,-1.9) 


31.4 


14.2 ± 0.2 ± 1.0 


(+0.8,-0.7) 


37.5 


4.7 ±0.1 ±0.3 


(+0.2,-0.3) 


43.6 


1.64 ±0.06 ±0.09 


(+0.10,-0.12) 


50.0 


0.59 ±0.03 ±0.03 


(+0.03,-0.04) 


61.3 


0.120 ±0.010 ±0.008 


(+0.009,-0.009) 


82.7 


0.0109 ±0.0024 ±0.0009 


(+0.0011,-0.0013) 



Table 1: Measured inclusive jet cross-section da/dEj, . The statistical and 
systematic uncertainties —not associated with the absolute energy scale of the jets- 
are also indicated. The systematic uncertainties associated to the absolute energy 
scale of the jets are quoted separately. The overall normalization uncertainty of 
2.25% is not included. 
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M 


(E^) 4 (E iet d 3 a/d^d^Y d ^z)v ± stat - ± s y st - 


syst. E'^-scale 




(xO.389 10 6 ) 


(xO.389 10 6 ) 



{W 1P ) = 180 GeV 



0.209 


18.7 ±0.2 ±1.1 


(+0.7,-0.7) 


0.246 


12.3 ±0.2 ±0.4 


(+0.5,-0.6) 


0.288 


7.38 ±0.15 ±0.49 


(+0.36,-0.37) 


0.346 


4.14 ±0.15 ±0.45 


(+0.23,-0.24) 


0.465 


1.36 ±0.06 ±0.17 


(+0.07,-0.09) 



(W 7P ) = 255 GeV 



0.209 


10.2 ±0.2 ±0.4 


(+0.4,-0.4) 


0.246 


6.86 ±0.20 ±0.40 


(+0.36,-0.44) 


0.288 


4.41 ±0.17 ±0.22 


(+0.12,-0.20) 


0.346 


2.12 ±0.15 ±0.14 


(+0.19,-0.10) 


0.465 


0.66 ±0.06 ±0.08 


(+0.04,-0.05) 



M 


Ratio ± stat. ± syst. 


syst. -scale 


0.209 


1.21 ±0.03 ±0.09 


(+0.01,-0.00) 


0.246 


1.19 ±0.04 ±0.07 


(+0.01,-0.02) 


0.288 


1.11 ±0.05 ±0.08 


(+0.02,-0.00) 


0.346 


1.30 ±0.10 ±0.11 


(+0.04,-0.01) 


0.465 


1.37 ±0.13 ±0.11 


(+0.01,-0.02) 



Table 2: Measured scaled jet invariant cross section 

(E^ t ) 4 {E^ et d 3 a/dp l ^dp l Y t dp l ^ t ) TI as a function of iy for (W yp ) = 180 and 255 
GeV and their ratio after correcting for the difference in the photon flux. Other 
details as in the caption to Table 0. 
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El£ range 
(GeV) 


a s (M z ) ± stat. 


± syst 


. ± th. 


<^ Ct > 

(GeV) 


a s ((£'^ !t )) ± stat. 


± syst. ± th. 


17 


- 21 


0.1229 ±0.0001 


+0.0026 
-0.0017 


+0.0055 
-0.0042 


18.6 


0.1644 ±0.0003 


i n nn/i c 
+U.UU40 

-0.0031 


+U.U1UO 
-0.0076 


21 


-25 


0.1217 ±0.0002 


+0.0020 
-0.0023 


+0.0052 
-0.0043 


22.7 


0.1558 ±0.0004 


+0.0033 
-0.0038 


+0.0087 
-0.0071 


25 


-29 


0.1210 ± 0.0004 


+0.0026 
-0.0032 


+0.0053 
-0.0041 


26.7 


0.1497 ± 0.0006 


+0.0040 
-0.0049 


+0.0083 
-0.0064 


29 


-35 


0.1228 ±0.0006 


+0.0030 
-0.0039 


+0.0058 
-0.0048 


31.4 


0.1477 ±0.0008 


+0.0044 
-0.0057 


+0.0086 
-0.0069 


35 


-41 


0.1228 ±0.0011 


+0.0040 
-0.0034 


+0.0069 
-0.0053 


37.5 


0.1429 ±0.0015 


+0.0056 
-0.0046 


+0.0095 
-0.0072 


41 


-47 


0.1194 ±0.0025 


+0.0058 
-0.0052 


+0.0081 
-0.0064 


43.6 


0.1347 ±0.0032 


+0.0074 
-0.0066 


+0.0105 
-0.0082 


47 


-55 


0.1227 ±0.0043 


+0.0064 
-0.0061 


+0.0105 
-0.0074 


50.0 


0.1355 ±0.0053 


+0.0079 
-0.0074 


+0.0130 
-0.0091 


55 


- 71 


0.1186 ±0.0098 


+0.0092 
-0.0111 


+0.0134 
-0.0093 


61.3 


0.1262 ±0.0112 


+0.0106 
-0.0126 


+0.0153 
-0.0106 



Table 3: The a s (Mz) values determined from the QCD fit of the measured 
da/dE^ in the different E^ regions and the a s ((E } f: t )) values determined as a 
function of E^ . The statistical, systematic and theoretical uncertainties are also 
indicated. 



17 



ZEUS 



ZEUS 98-00 
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LO QCD 



142 < W < 243 GeV 
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Figure 1: a) Measured inclusive jet cross section, da/dE 3 ^ (dots). The thick 
error bars represent the statistical uncertainties of the data, and the thin bars show 
the statistical and systematic uncertainties added in quadrature. The uncertainty 
associated with the absolute energy scale of the jets is shown separately as a shaded 
band. The LO (dashed line) and NLO (solid line) QCD parton-level calculations 
corrected for hadronisation effects are also shown, b) The fractional difference 
between the measured da/dE^ and the NLO QCD calculation; the hatched band 
shows the uncertainty of the calculation. 
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ZEUS 




Figure 2: Measured scaled jet invariant cross section, 

(E 1 ^ t ) 4 (E^ et d 3 a/dp l xdpY t dp i ^ t } v , averaged over —2 < rf^ < 0, as a function 
of x T (dots) for a) (W 1P ) = 180 GeV, b) (W 1P ) = 255 GeV. Other details are as 
given in the caption to Fig. ^ 
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ZEUS 



• ZEUS 98-00 

NLOQCD 

scaling 



<W > = 180 GeV/<W > = 255 GeV and -2 < i/*< 

TP TP 'TP 




jet energy scale uncertainty 
EZZZZ] NLO uncertainty 



0.2 



0.3 



0.4 



0.5 



Figure 3: Measured ratio of scaled jet invariant cross sections, after correcting 
for the difference in the photon flux between the two W 1V intervals, as a function 
of xt (dots). The dashed line is the scaling expectation. Other details are as given 
in the caption to Fig. El. 
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Figure 4: a) The a s (Mz) values determined from the Q CD fit of the measured 
da/dE 3 ^ 1 in the different E 3 ^ regions (open circles). The combined value of a s (Mz) 
obtained using all the E 3 ^ regions is shown as a dot. b) The a s (E^ t ) values de- 
termined from the QCD fit of the measured dajdE 3 ^ as a function of E 3 ^ (open 
circles). The solid line represents the prediction of the renormalisation group equa- 
tion obtained from the a s [Mz) central value as determined in this analysis; the 
light-shaded area displays its uncertainty, c) The l/a s (E 3 ^ ) values as a function 
of E^ (open circles). The solid line represents the result of the two-loop a s fit 
to the measured values. The dashed line represents the extrapolation of the result 
of the fit to E 3 ^" = Mz- The dot, plotted at E 3 £ = Mz, represents the inverse 
of the combined value shown in a). In all figures, the inner error bars represent 
the statistical uncertainties of the data and the outer error bars show the statistical 
and systematic uncertainties added in quadrature. The dashed error bars represent 
the theoretical uncertainties. The current world average [32] (dotted line) and its 
uncertainty (shaded band) are displayed. 
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